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Abstract 

QCD corrections to the electroweak one-loop result for the partial 
width T(H — > 66) are studied. For the decay channel into bottom 
quarks the rate is affected by a virtual top quark through electroweak 
interactions. The calculation of QCD corrections to this quantity is 
performed for an intermediate range Higgs mass in the heavy top 
mass limit. The leading correction of order 0(GFa s mf) is estimated. 
Numerically the contribution is of comparable size as the electroweak 
correction, but of opposite sign. 
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1 Introduction 



At present high energy colliders experimental data have turned out to be in re- 
markable agreement with theoretical predictions of the Standard Model. Precision 
experiments have covered many aspects of both its electroweak and QCD sector. 
Although the formulation of the Standard Model as a SU(3)c x SU(2)l x U(1)y 
gauge theory is firmly established, experimental evidence for the physical Higgs 
boson is still missing. Since the Higgs mechanism of spontaneous symmetry break- 
ing and mass generation would reveal its structure through the properties of the 
Higgs boson, the high mass Mh of this scalar particle hides the Higgs sector com- 
pletely from observation at present energies. Future accelerators like LHC and 
NLC, which are dedicated to the search and the study of the Higgs particle, may 
hopefully close the energy gap to its production threshold. 

For intermediate Higgs masses Mjj < 2M\y the dominant decay channel is the 
decay into bottom quarks. A theoretical prediction for the partial width T(H — > 
bb) is therefore needed for the physics analysis with high precision. Consequently 
much work has been spent in the past on the calculation of radiative corrections 
to this quantity. Excellent reviews on the Higgs phenomenology can be found for 
example in [2J. 

Electroweak radiative corrections, which were studied at the one-loop level by 
several groups |3], |J are of particular interest for the process H — > bb, since 
they are affected by the mass of the top quark. Virtual top states are possible due 
to bottom-top transitions mediated by charged Higgs ghosts ^ or W ± bosons. 
The leading m 2 dependence originates from diagrams with an exchange of $ ± 
due to the Yukawa coupling of the Higgs ghosts to the fermion line. 

In this work QCD corrections to the electroweak one-loop result are calculated 
in the heavy top mass limit. For this purpose the hard mass procedure [§, |7|, |[ || 
is employed. It results in a power series in the inverse heavy mass l/m t , of which 
we compute the leading term of order 0(GFa s m^). 

It is convenient to consider the corresponding scalar current correlator n(g 2 ), 
since its absorptive part completely determines the partial decay rate 

T(H^bb) = ^-lmIl(Mj I ). (1) 

The calculation of n(g 2 ) involves multiloop propagator type integrals as well as 
massive tadpole integrals. Many of the computational tools, mostly based on the 
algebraic manipulation language FORM [|1(J, have already been used previously 
for the calculation of 0[Gpa s vn? t ) corrections to the partial width T(Z — > bb) of 
the Z boson JIT] . Nevertheless, the structures of the calculations for the processes 
H ^ bb and Z — > bb differ in various respects. First, the generic diagram with 
a primordial decay of the Higgs into a pair of charged Higgs ghosts and their 
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Figure V. The electroweak diagrams contributing to the self energy of the H 
boson. Internal dashed line: Higgs ghost, thin lines: bottom quark, thick lines: 
top quark. 

subsequent transition into b quarks does not contribute to the considered order. 
This can be seen from the H& + &~ coupling and is explained on dimensional 
grounds. Second, the Yukawa couplings of the Higgs to the scalar fermion cur- 
rents are proportional to the quark masses. Combined with mass terms from the 
fermion traces the decay rate is proportional to the square of the bottom mass 
and would vanish in the limit m& = 0. Third, the quadratic dependence of the 
correction term in the considered order on the bottom as well as the top mass 
affects the renormalization of the scalar correlator. Wheras for the Z decay the 
sum of the three loop diagrams was finite, this is not the case for the Higgs rate. 
The renormalization of the quark masses induces contributions from lower order 
diagrams which have to be included to arrive at a finite result. 

In Section 2.1 the calculation of the three loop propagator diagrams is de- 
scribed. Subsequently the relation between renormalized and bare quark masses 
including corrections up to the two loop level is combined with the Born graph 
as well as pure electroweak and QCD diagrams. The corresponding contributions 
to the order 0(GFa s mf) corrections are discussed in Section 2.2. In Section 3 
universal corrections to T(H — > bb) are added to the flavour specific corrections. 
The numerical size of the effects is evaluated. 

2 Description of the Calculation 

2.1 Three-Loop Diagrams 

All diagrams of order 0(GF<y s mf) are obtained from two different generic graphs 
(see Figure 1), to which a gluon has to be attached in all possible ways. The 
calculation is performed with dimensional regular isat ion in the on-shell scheme. 
In D = 4 — 2e dimensions the definition of the Hermitian, anticommuting 75 
matrix is used. After combining pairs of 75 to unity, traces with single 75 vanish 
identically. We work in the t'Hooft gauge with an arbitary QCD gauge parameter 
£s. The verification of gauge invariance serves as an internal check of our result. 

Being interested in the leading correction of the power series with respect to 
the inverse heavy mass l/m t , we apply the hard mass procedure as developed 
and elaborated in || [?J |8|, ||. Those subgraphs are selected, which comprise all 
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top propagators and become one-particle-irreducible after contracting the heavy 
lines to a point. They are expanded with respect to external momenta and small 
masses. This formal Taylor expansion is then inserted as an effective vertex into 
the remaining diagram which has to be evaluated. In the heavy top limit the 
bottom mass as well as the W mass and the Higgs mass are considered small 
as compared to the top mass. Since the expansion is performed only to the first 
nonvanishing order, the Higgs ghost propagators are always contained in the hard 
subgraphs. For the leading term My/ is therefore set to zero. As a consequence 
the electroweak gauge parameter drops out trivially. Although the bottom mass 
is also a small expansion parameter, it cannot be neglected throughout, since that 
would result in an identically vanishing decay rate. The bottom mass is not only 
introduced as an overall factor through the Yukawa coupling, but is also present 
in fermion traces with an odd number of Dirac matrices for mj = 0. Besides the 
expansion of the hard subgraphs as prescribed by the hard mass procedure also 
the remaining diagram is expanded with respect to m\. After all expansions are 
performed, the leading piece can be isolated. Its contribution AT aas to the partial 
Higgs decay rate reads 

Ar QQs = T m 2 b x t ^ j<4 + - f 19 + 12 In 4 - 15 In ^f] 

199 . . 2 3 o , u 2 M 2 H , s 

+ — -24C 3 )-— 7r 2 + 161nAr -411n— £ ( 2 ) 
3 6 mf /r 

+9 In At H In 2 —S. - 18 In At In ■ 1 



with 



m 2 2 fi 2 mf /i 2 



( "" N C M H 



4\/27r 
Gftti 2 



(3) 



8\/27r 2 ' 

According to the conventions of the multiloop integration package MINCER Jl2[ , 
which we used for the calculation of massless integrals, terms with 7^ and ln(47r) 
are suppressed. In addition for each loop integration a factor e^ 2 ^ 2 / 2 is included 
for convenience. It is clear, that physical results are not influenced by this con- 
vention, since only terms of order e 2 are affected. A consistent convention holds 
also for the calculation of massive tadpole integrals [|l]| . Our result is obtained in 



the limit of large top masses and is applicable only in the range of intermediate 
Higgs masses, in particular in the regime Mu < 2m t below the top threshold. 
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2.2 Mass Renormalization and Induced Contributions 



mY = mi, 



The result of eq.(|2|) is still divergent. Renormalization of the quark masses induces 
additional contributions which lead to a finite expression for the Higgs decay rate. 
The relation between the bare mass m b and the renormalized mass m° s of the 
bottom quark in the on-shell scheme may be written as follows 

1 + C a x t + C as ^ + C aas x t ^\ (4) 

7T 7T J 

where only corrections of orders relevant for our problem are taken into account. 
The corresponding coefficients C a , C as and C aas can be derived from the fact that 
the renormalized (pole) mass m° s is defined through the location of the pole of 
the quark propagator 

Sf «m k -|l + S(j))' 5 
The selfenergy of the bottom quark is decomposed in the form 



S(p) = m b S 2 [ 1 + (m b - iS)T, x . (6) 



The functions Si^ receive contributions of the order 0(a s ) from diagram 2a 
with 5 = 1, of order (9(GVm 2 ) from diagram 2b with 5 = (1 — 75) and of order 
0(GF<y s mf) from diagrams 2c-2e with 5 = (1 — 75). For the calculation of the 
quark selfenergies we again employ the hard mass procedure. As suggested in fT4| 
remaining standard scalar integrals may be simplified by an expansion in m 2 /p 2 
around 1. Neglecting all higher order terms only integrals on the bare mass shell 
need to be evaluated: 

m° s = m b [1 + E 2 (l) - E 2 (l) (2^(1) + E x (l)) + • ■ ■] (7) 

The derivatives S 2 = <9£ 2 /<9(m 2 /p 2 ) ma Y be conveniently obtained through 
derivations with respect to the bare mass, thus raising the power in the denomi- 
nator of the quark propagator. One obtains 

2 
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Figure 2: Diagrams for the renormalization of the bottom quark. 



A replacement of the bare bottom and top masses through the renormalized ones 
according to eq.([|) in the results for the Born graph, the electroweak and the 
QCD corrected diagrams leads to induced contributions to the order 0(Gjra; s mf). 
Combined with eq.(|]) one obtains the finite result for the decay rate. 



3 Discussion 



In the previous section the non- universal corrections to the partial width T(H — > 
bb) have been discussed. Combining all contributions one obtains the following 
result for these process dependent corrections: 



a -pnon—univ. 
ZAi H^bb 



-6x t T ml 



1 + 



O.S 



2 In 



M 



H 



7T 



mi 



(9) 



This result has been confirmed recently in JTSI . 

In addition the partial rate is also influenced by universal corrections. Process 
independent terms of order OiGpOigm 2 ) originate from Ar due to the use of the 
Fermi constant Gp instead of the electroweak coupling constant a and from the 
renormalization of the Higgs vacuum. They are given by |4| 



a -puniv._ 
lAL H^bb 



-T m 2 b 



n 



WW 



(0) 



+ Ren™(M£) + ■ ■ 



M 2 W 



1 , r a s 
1 + OqcD 



7T 



(10) 



where the dots indicate terms of subleading order and the QCD correction factorQ 
is given by 5 QC d = 3-2 ln(M|/mg). 

For the calculation of the the (unrenormalized) selfenergies n WM/ and U HH 
again the hard mass procedure can be applied. Being interested in their real parts 
the leading hard subgraphs are constituted by the two loop diagrams themselves 
with the W selfenergy to be evaluated for zero external momentum. Since the 
hard mass procedure for the Higgs vacuum polarization is in fact an expansion 
of the graph with respect to the external momentum, the derivative Rell^^' is 
readily obtained from the next to leading contribution to the power expansion. 



1 We thank the authors of for pointing out that in the earlier version of this paper 
this factor was missing in eq.(nQ). As a consequence the interference term of the 1-loop QCD 
corrections and the 1-loop universal electroweak corrections was erraneously not taken into 
account. 
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After top mass renormalization the combination eq.fllOD of the W- and the H 
selfenergy is finite in their sum at order 0(GFQ. s m 2 ). The result reads 



a -pumv._ 
^ L H^bb 



T m 2 b x t 



7-2- 



a. 



71 



1 + 5, 



a. 



QCD' 



71 



The electroweak term in the first bracket is in agreement with |4], |5j and its QCD 
contribution reproduces the result obtained by Jl5| (see also |16|, [17|). 
In the sum we arrive at the final result for the partial decay rate 



H^bb 



r ml I 1 + x t 




-1 - -7l 2 



01 Ml 

2 log H 



m 2 b 



(12) 



The QCD correction to the electroweak result is of comparable size as the elec- 
troweak correction itself, but of opposite sign. With m b = 4.7 GeV, m t = 174 
GeV and M H = 60 GeV (120 GeV) the electroweak contribution of 0.32% nor- 
malized to the Born width is combined with a QCD corrected term of —0.21% 
(-0.24%). 

To conclude, in this work we have calculated the QCD corrections to the 
electroweak one-loop result for the partial Higgs decay rate T(H — > bb) in the 
heavy top mass limit. For intermediate Higgs masses electroweak corrections are 
significantly reduced by the 0(Gpa s m 2 ) corrections. 
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